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ABSTRACT 

We have serendipitously identified the first lithium-rich giant star located close to the red giant 
branch bump in a globular cluster. Through intermediate-resolution FLAMES spectra we derived a 
lithium abundance of A(Li)=2.55 (assuming local thermodynamical equilibrium), which is extremely 
high considering the star’s evolutionary stage. Kinematic and photometric analysis confirm the object 
as a member of the globular cluster NGC 362. This is the fourth Li-rich giant discovered in a globular 
cluster but the only one known to exist at a luminosity close to the bump magnitude. The three 
previous detections are clearly more evolved, located close to, or beyond the tip of their red giant 
branch. Our observations are able to discard the accretion of planets/brown dwarfs, as well as an 
enhanced mass-loss mechanism as a formation channel for this rare object. Whilst the star sits just 
above the cluster bump luminosity, its temperature places it towards the blue side of the giant branch 
in the colour-magnitude diagram. We require further dedicated observations to unambiguously identify 
the star as a red giant: we are currently unable to confirm whether Li production has occurred at the 
bump of the luminosity function or if the star is on the pre zero-age horizontal branch. The latter 
scenario provides the opportunity for the star to have synthesised Li rapidly during the core helium 
flash or gradually during its red giant branch ascent via some extra mixing process. 

Subject headings: globular clusters: individual (NGC 362) - stars: abundances - stars: Population ii 


1. INTRODUCTION 

Lithium is one of the most fragile elements, promptly 
consumed by (p,a) reactions when the temperature in 
stellar interiors reaches Ta:2.5xl0 1 2 * 4 5 6 K. According to 
standard stellar evolution theory, when a star leaves 
t he main-se quence and experiences the first dredge-up 
ijlbenl Il967h . the convective envelope penetrates into 
hot internal layers, mixing the outer zones with mat¬ 
ter that has undergone partial H-burning. The main 
surface abundance changes are increases in 4 8 * He, 14 N 
and 13 C , and a decrease in the 12 C abundance (e.g., 
iKarakasI I2010H . The initial Li content is reduced by 
factors 15—20 in Population II stars, reaching typical 
values of A(Li) <1.5 dex. Observations of Li abun¬ 
dances in low-mass red giant-branch (RGB) stars have 
indeed provided strong and unambigu ous con fir mat ion to 
these theoretical pr edictions (e.g.. lLambert et ahlll98Ct 
IGratton et a. 1.1 [200011 . Yet, approximately between 1—2 
% of K giants have been found to exhibit significantly 
higher Li abundances and they are thus referred to as 
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Li-rich giants (1 Brown et al.l 1 19891 ^ Monaco et all 120111 : 


iKumar et all [201 lb iLebzelter et all 20121 and references 
therein). Interestingly, this fraction might be even 
l ower, that is 0.56%±0.16%, as recently emphasised by 
IZhang et all (l2015i ). A clear explanation for this rare 
phenomenon is still wanting and different mechanisms 
have been suggested. They envisage co n tam ination from 
ejecta of nearby novae (Martin et all Il994fl. ac cretion 
of planets and/or brown dwarfs (jSiess fc Livfol Il999(h 
as well as freshl y synthesised Li in the star s them¬ 
selves (see e.g., iSackmann fe Boothrovdl I1999T) . The 
physical pr ocess responsible fo r Li p roduction is the 
well known iCameron &; Fowleil (j 197 lh mechanism: the 
fusion sequence is 3 He(a, 7 ) 7 Be(e,i') 7 Li, with the re¬ 
quirement that the astrated 7 Be is quickly trans¬ 
ported to the cooler envelope so that the eventual 7 Li 
produced will n ot be destr oyed by proton captures. 
iCharbonnel fe Balachandraiil (j2000f ) proposed two dis¬ 
tinct episodes of Li production depending upon the 
mass of the star. In low-mass RGB stars, with de¬ 
generate helium cores, Li production may be associated 
with the bump in the luminosity function; alternatively, 
intermediate-mass stars can synthesise Li once the core- 
He burning has ceased and they ascend the early asymp¬ 
totic giant branch (AGB, see that paper for details). 

The preliminary observational detections of Li-rich gi¬ 
ants seemed to suggest that these stars occupy only spe¬ 
cific regions of the Hertzsprung-Russel diagram (such as 
the RGB bump, tip and/or the red clump). However, 
sever al studies (e.g-. lMonaco et aJ1l2011h ILebzelter et al.1 
f201S| ) have recently identified Li-rich stars at magnitudes 
that span the red giant branch. Their results indicate 
that the formation of Li-rich giants might not be limited 
to the aforementioned specific evolutionary events. 

Given how rare these stars are, large surveys are 
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required (and underway) to increase the number of 
observed Li-rich objects. More detections are nec¬ 
essary to improve our understanding of the physi¬ 
cal phenomena involved and provide empirical con¬ 
straints on their formation. Many stellar environments, 
within and beyond the Milky Way, have been scru¬ 
tinised for the presence of Li-rich giants (thin/thick 
disk, bulge and halo in the Milky Way, dwarf ga lax¬ 
ies, o p en and globular cluste r s -GCs, e.g.. iMonaco et al.1 
20111 iLebzelter et all 120121 iMart ell fc Shetronel 2 0131 


Kirby et al. 2012i : iMonaco et alJl2014l : iPilachowski et ahl 

20001) . Clusters (either open or globular) offer an obvi¬ 
ous advantage in that they host mono-metallic, almost 
coeval stellar populations and hence represent an excel¬ 
lent target for this kind of investigation: their masses 
and luminosities are also much more accurately deter¬ 
mined with respect to field stars. To date, only three 
Li-rich giants have been reported in GCs: V42 in M 5 
dCarnev et al.lll998D. V2 in NGC 362 (iSmith et al .1119991) . 
and IV-101 in M3 feraft et al.lll999Tl . All of them are lo¬ 
cated in the upper part of the RGB: V2 and IV-101 lie 
at the tip of the RGB, while V42 is actually a W Vir 
post-AGB star. We present in this Letter the serendip¬ 
itous discovery of a Li-rich giant close the RGB bump 
in the globular cluster NGC 362, a relatively metal-rich 
cluster ([Fe/H]=—1.26 dex, 2011 update of the Harris’ 
catalogue), characterised b v chaotic and unusual orbital 
parameters (|Dinescu et al.lfl999D . This is the first time 
that a Li-rich giant is detected at this evolutionary stage 
in a GC. 


2. OBSERVATIONS AND ANALYSIS 

Observat ions were conducted in visitor mode with 
FLAMES (|Pasauini et al.ll2002j ) mounted at ESO VLT 
on 11, 12, and 13 December 2014, under program 094.D- 
0276(A). Spectra were collected with the scientific aim 
of inferring Li and Al abundances for large samples 
of RGB stars in sever a l GCs, as pa rt of our dedi¬ 
cated survey (see iD’Orazi et al.l 12014 hereinafter Pa¬ 
per i). We utilised the HR15N setup that provides a 
nominal resolution of R=17,000 and a spectral cover¬ 
age from 6470 A to 6790 A, allowing the simultane¬ 
ous inclusion of the strong doublets for Li I (6708 A) 
and Al I (6696/6698 A). The data were reduced using 
the dedicated ESO software (version 2.12.2, available at 
http://www.eso.org/sci/software/pipelines/) and 
resulting in bias-subtracted, flat-fielded, wavelength- 
calibrated, one-dimensional spectra. Sky subtraction 
and rest-frame shift were then carried out within IRAlfl 
As with Paper I, stellar parameters were derived in the 
following manner: effective temperature s (7## are in¬ 
ferred fr om V and K, photom etry ffrom iMomanv et all 
12004 and iSkrutski e et al.ll 2006l re spectively') and the cal¬ 
ibration bv lAlonso et all (I1999D . The reddeni n g and 
metallicity values are adopted from the lHarrisI (1996) 
catalogue, that is E (B — V)=0.05 and [Fe/H]=—1.26 
dex. Surface gravities were computed from the stan¬ 
dard formula relating luminosities, temperatures and 



Wavelength (A) 


Fig. 1.— Spectral synthesis of the Li I line at 6707.78 A for our 
target star 

massct!3, adopting M=0.85 M 0 and bolometric solar 
magnitude of M{, o ; )0 = 4.75. Microturb ulence veloci¬ 
ties (£) are from the relation bv iGratton et al.l ( 1996 ). 
Abundances were obtained via spectral synthesis for both 
Li and Al using the local thermodynamical equilibrium 
(LTE) code MOOG (lSnedenlll973l . 2014 version); interpo¬ 
lated Kurucz model atmospheres based on the ATLAS9 
cod e with q-enhancement an d no convective overshoot¬ 
ing (jGastelli fc Kuruc^l2004 ) were used throughout this 

studjO 

We adopted the same line list as in Paper I, perform¬ 
ing the comparison between synthetic and observed spec¬ 
tra between 6696 and 6720 A (exploiting the Ca I line 
to evaluate the spectral broadening). Errors related to 
best fit determination and uncertainties in stellar param¬ 
eters were consistently calculated as for Paper I and then 
added in quadrature, given their independence. The to¬ 
tal internal random error results in 0.09 dex and 0.12 dex 
for Li and Al, respectively. Finally, non-LTE corrections 
we re applied to Li abundances, following prescriptions 
bv lLind et al.l (120091) . 

3. RESULTS AND DISCUSSION 

Stellar parameters, kinematics and abundances for Li 
and Al are reported in Table [T] for our target star. We 
inferred a Li abundance of A(Li)=2.55±0.09 (0.13 dex 
lower when departures from LTE are taken into account) 
which is significantly higher than our determined aver¬ 
age for the cluster (A(Li)=0.90±0.02 rms=0.15, 66 RGB 
stars; D’Orazi et al. 2015, in preparation). Thus, star 
#15370 joins the small family of Li-rich giants in GCs 
whose other members include a star in M3, M5 as well as 
a neighbour also in NGC 362. Star #15370 distinguishes 
itself because it is the first Li-rich giant discovered in a 
GC near the luminosity of the RGB bump (roughly 0.12 


9 IRAF is the Image Reduction and Analysis Facility, a general 
purpose software system for the reduction and analysis of astro¬ 
nomical data. IRAF is written and supported by National Optical 
Astronomy Observatories (NOAO) in Tucson, Arizona. 
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magnitudes brighter than Vh Mmp =15.40 from iNataf et al.l 
120131 ). Its siblings all sit close to the tip of their giant 
branches. There is little doubt that the star is Li rich: the 
intrinsic strength of the 6708 A doublet is evident from 
Figure [1] where we show the observed spectrum along 
with synthetic profiles calculated for different Li abun¬ 
dances. The equivalent width (EW), measured using the 
IRAF task splot and a gaussian profile, is EW=250 mA. 
Unfortunately, our spectral coverage does not include the 
other Li I feature at 6104 A, thus we cannot utilise this 
line as a secondary indicator. 

Kinematic information confirms cluster membership. 
The star’s radial velocity V ra d= 222.68±1.02 km s -1 is 
consistent with the cluster average we derived from our 
sample of 71 RGB stars, i.e., \ r ad= 222.78±0.62 km s -1 . 
This value a grees ve r y w ell with V ra[ j=223.5±0.5 km s _1 
as given bv IHarrisl (119961) . Moreover, the position on 
the colour-magnitude diagram (Figure [2]) as well as the 
spectral comparison with another cluster member with 
almost identical stellar parameters (see Figure 3: namely 
star #2213, T eff =4821K, logg=2.19 cm s" 2 , £=1.51 
km s -1 ) further advocate that this Li-rich giant belongs 
to NGC 362. 

We find an A1 abundance of [A l/Fe1=0.08 dex, which 
according to the definition by e.g.. lCarretta et all (|2009t l 
and Paper I suggests that our Li-rich giant is a member of 
the first generation of t he cluster po pulatio rPl . A sim¬ 
ilar result was found bv lSmith et al.l ( 1999 ) who derived 
high O and low Na and A1 abundances for their Li-rich 
giant V2. Both IV-101 and V42 seem to show a (mod¬ 
est) Na enrichment and probably belong to the second 
generation of their respective clusters. Given the quite 
limited sample we are currently dealing with, statistics 
on the occurrence of Li-rich stars among the different 
stellar populations in GCs is not meaningful. We hope 
that the dedicated surveys will lead to a higher number 
of detections and help with such analyses. 

One possible explanation for the formation of Li-rich 
K giants is by accre tion of planets and/or brown dwarfs 
(iSiess fo Livio!fl999 ). Given the high Li abundance of 
A(Li l)«2.5 dex, this scenario seems unlikely for our star 
because it would require the ingestion of a very massive 
companion. Giants during this phase of evolution typi¬ 
cally have a convective envelope of ~0.5 M 0 . To account 
for the observed Li abundance, star #15370 needs to have 
engulfed a companion wit h mass o f at least rH.O M 0 (us¬ 
ing Equation 2 given bv lSiess fc Livioll 19991 and adopting 
a Li content for the donor star of 2.55 dex, which is al¬ 
ready an unusually high value with respect to the Spite 
plateau). Such a big accretion event should also result 
in increased rotational velocity and /or e nha nced level 
of chromospheric activity (see e.g., iHurlev et al.l 120021 : 
iDrake et al.ll2002l ) , which is not evident from our spectra 
(at least within the limits imposed by the quality of our 
da ta). _ __ 

Ide La Reza et alJ (U996U1997H proposed a model con¬ 
necting the high Li abundances observed in K giants to 
the evolution of circumstellar shells. In this framework, 
the internal mechanism responsible for the Li enrichment 

12 A discussion on the multiple population phenomenon in GCs 
is not the purpose of the current manu script; we refer the reader 
to recent reviews by IGratton et alJ l|2012l 'l an d IPiotto et alJ (|2012H 
for spectroscopic and photometric perspectives, respectively. 


will initiate a prompt mass-loss event. Empirical evi¬ 
dence for this paradig m has not been forthcoming (e.g., 
iFekel fe Watsonlll998l ). We inspected our spectra around 
the Ha feature, looking for possible blue-shifted asym¬ 
metry as a hint of enhanced (gas) mass-loss. Within the 
limit of our S/N and resolution, we did not discover any 
asymmetry in the Ha profile of star #15370 compared 
to other Li-normal giants with similar atmospheric pa- 
ra meters. In this respe ct, w e confirm results publ ished 
bv lMonaco et al.l (1201 11 1 and iLebzelter et all (|2012f h 

The literature offers two alternative scenarios that are 
consistent with our object: 

(i) Star #15370 belon gs to the class of Li-rich 
giants described by ICharbonnel fe Balachandranl 
(200 0) that are synthesing Li at the beginning 
of the RGB bump phase. After the hydrogen¬ 
burning shell removes the composition disconti¬ 
nuity left behind by first dredge-up, some extra¬ 
mixing process facilitates Li production via the 
Cameron & Fowler mechanism. Subsequently, 
once the mixing deepens enough to convert ad¬ 
ditional 12 C to 13 C, the surface material is ex¬ 
posed to temperatures that are too high and the 
fresh Li is quickly destroyed. Thus, the Li- 
rich phase is swift and a star will not retain its 
peak Li abundance o nce the 12 C/ 13 C dr ops belo w 
the s tandard value (jCharbonnel fe Balachandranl 
12000 1. Parametrised analyses, unsurprisingly, 
require that these stars have unusually large 
diffusion coefficie nts to allow f or the rap id 
transp ort of ‘ B e dSackmann fe Boothrovdl 119991 : 
iDenissenkov fe Herwigll2004l l. 

(ii) A second possibility is that the production of 
Li mig ht hav e occu rred during the helium flash 
dKumar et al.l l20lHl or gradually as the star ap¬ 
proached the RGB trp_ d ue t o some kind of extra 
mixing (Latt anzio et alJ 2Q15|) . Asteroseismic anal - 
ysis of KIC 5000307 by Silva Aguirre et all (120141 1 
has identified this Li-rich field star as clearly un¬ 
dergoing core helium burning. An analogous sce¬ 
nario for star #15370 would require that it is ac¬ 
tually a pre zero age horizontal branch (ZAHB) 
star. A close inspection of the CMD indicates that 
the star is located on the blue side of the RGB 
and its position is not inconsistent with theoret¬ 
ical predict i ons o f pre-ZAHB evolution published 
bv ISilva Aguirre et all (120081 ). 

As a further check we have run evolutionary models 
using MONSTAR (Monash Versio n of the Mt. Stromlo 
evolution code, see lAngelou et al.ll2012lh For our model 
we assumed a metallicity of [Fe/H]=—1.26 dex, an a- 
element enhancement of [a/Fe]=0.2, and initial mass of 
0.81 M 0 corr espon din g to a main-sequence turn off age 
of 11.5 Gyr (iDotter et al.l 1201 71) . As can be see from 
Figure [5] the star’s location is best explained by RGB 
evolutionary tracks. We note that although the star is 
(slightly) on the blue side of the giant branch it is still 
somewhat removed from the cluster HB. However, in or¬ 
der to unequivocally determine the evolutionary phase, 
we plan to acquire high-resolution (R«40,000), high S/N 
spectra (>100) of our target star along with one of its 
siblings to derive their carbon isotopic ratios 12 C/ 13 C. In 
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Fig. 2.— Colour-magnitude diagram for NGC 362 (Momany et 
al. 2004). The target star is emboldened. 


principle, differences in Fe II lines would also provide us 
with evidence of a slightly different gravity for $T5370 
an indication that the star experienced mass loss at the 
tip of the RGB branch (a difference in mass of approx¬ 
imately 0.2 Mq implies a difference in gravity of about 
0.10 dex). 
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Fig. 3.— Spectral comparison for stars #15370 and #2213. 


TABLE 1 

Stellar parameters, kinematics and chemical information for star #15370 


Star ID 

RA DEC 

(hh:mm:ss) 

B 

(mag) 

V 

(mag) 

K s 

(mag) 

V rad 

(km s -1 ) 

S/N 

(@6708A) 

T e ff 

(K) 

log g 
(cm s~ 2 ) 

e 

(km s -1 ) 

A(Li)LTE 

dex 

A(Li)NLTE 

dex 

[Al/Fe] 

dex 

15370 

01:03:43.9 -70:52:20.9 

16.052 

15.233 

12.944 

233.47 

190 

4812 

2.18 

1.52 

2.55 

2.42 

0.08 















